Adsorption of H2O, NH3, CO, NO2, and NO on graphene: A first-principles study 
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Motivated by the recent realization of graphene sensors to detect individual gas molecules, we 
investigate the adsorption of H2O, NH3, CO, NO2, and NO on a graphene substrate using first 
principles calculations. The optimal adsorption position and orientation of these molecules on the 
graphene surface is determined and the adsorption energies are calculated. Molecular doping, i.e. 
charge transfer between the molecules and the graphene surface, is discussed in light of the density 
of states and the molecular orbit als of the adsorbates. The efficiency of doping of the different 
molecules is determined and the influence of their magnetic moment is discussed. 

PACS numbers: 68.43.-h, 73.20.Hb, 68. 43. Be, 81.05.Uw 



I. INTRODUCTION 

The synthesis of monolayer graphite (i.e. graphene)^ 
and the experimental observation of Dirac charge carriers 
in this system^ ^ have awakened an enormous interest in 
this two-dimensional material. The unusual properties of 
carriers in graphene are a consequence of the gapless and 
approximately linear electron dispersion at the vicinity 
of the Fermi level at two inequivalent points of the Bril- 
louin zone. In the low-energy limit the quasiparticles in 
these systems are described in terms of massless chiral 
relativistic fermions governed by the Dirac equation. 

The good sensor properties of carbon nanotubes are 
already known for some time,^ but recently, the pos- 
sibility to use graphene as a highly sensitive gas sen- 
sor was also reported.^ It was shown that the increase 
in graphene charge carrier concentration induced by ad- 
sorbed gas molecules can be used to make highly-sensitive 
sensors, even with the possibility of detecting individual 
molecules. The sensor property is based on changes in 
the resistivity due to molecules adsorbed on the graphene 
sheet that act as donors or acceptors. The sensitivity of 
NH3, CO, and H2O up to 1 part per billion was demon- 
strated, and even the ultimate sensitivity of an individual 
molecule was suggested for NO2. These excellent sensor 
properties of graphene are due to two important facts: 
i) graphene is a two dimensional crystal with only a sur- 
face and no volume, which maximizes the effect of surface 
dopants, and ii) graphene is highly conductive and shows 
metallic conductance even in the limit of zero carrier den- 
sity. 

To fully exploit the possibilities of graphene sensors, 
it is important to understand the interaction between 
the graphene surface and the adsorbate molecules. We 
perform in this letter first principles calculations for the 
molecules NH3, NO2, NO, CO, and II2O adsorbed on 
graphene. We determine their exact orientation on the 
surface and their preferential binding site by calculat- 
ing their binding energy. Their charge transfer to the 
graphene surface is investigated in order to determine 
the donor or acceptor character of the molecular dopant. 



II. COMPUTATIONAL DETAILS 

The first principles calculations are performed using 
density funcional theory (DFT) which has been suc- 
cesfully used for the study of molecular a dsorbates on 
single-walled (carbon) nanotubes (SWNT).EEEEIIol All 
our DFT calculations were carried out with the Abinit 
code,Hl within the generalized gradient approximation 
(GGA) of Perdew, Burke, and Ernzerhof (PBE).^The 
advantage of GGA over the local density approximation 
(LDA) in this work is that the GGA will not lead to 
a strong bonding of the molecules as in LDA. So if the 
molecules bind in GGA, they will definitely bind in a real 
system (and in LDA) too. The distance between adsor- 
bate and the graphene surface, however, will be some- 
what overestimated and consequently the binding energy 
will be underestimated. 

We use a plane wave basis set with a cutoff energy of 816 
eV and pseudopotentials of the Troullier-Martins type^ 
For the adsorption of the molecules NH3, CO, and H2O 
we use non-spin-polarized calculations, while for NO2 
and NO, we use spin-polarized ones. The total system 
consists of a 4 X 4 graphene supercell (32 C atoms) with 
a single molecule adsorbed to it (Fig. [T]) and with a dis- 
tance of 16 A between adjacant graphene layers. The 
sampling of the Brillouin zone is done using a 5 x 5 x 1 
Monkhorst-Pack^^ grid. For the calculation of the den- 
sity of states (DOS) we use a 15 x 15 x 1 Monkhorst-Pack 
grid and a Gaussian smearing of 0.14 eV. 




FIG. 1: (Color online) H2O on graphene. 4x4 supercell of 
graphene with adsorbed H2O molecule. 

Charge transfers are calculated based on the Hirshfeld 
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charge analysis!^ The atomic charge Qa for each atom is 
obtained by (with p(r) the calculated density and p^iv) 
the electron density computed for the isolated atom A 
and taken from Ref. [TT]) 

from which the charge transfer (AQ) is deduced. From 
this result we determined whether or not the adsorbate 
acts as an acceptor or a donor. It should be noted that 
the size of the charge transfer is slightly dependent on 
the method used to calculate it. 

The distance from the adsorbate to the graphene surface 
is calculated from the difference in weighted averages of 
the different atoms of the molecule and the carbon atoms 
of the graphene sheet, where we used the atomic number 
Z of the atoms as the weight factor. 



III. RESULTS 

For each adsorbate three adsorption sites are consid- 
ered, namely on top of a carbon atom (T), the center of 
a carbon hexagon (C) and the center of a carbon-carbon 
bond (B) (see Fig. 1) . For these positions, different ori- 
entations of the molecules are examined and the adsorp- 
tion energy is calculated for all of them. The adsorption 
energy (E^) is the energy of the isolated graphene sheet 
and isolated molecule minus the energy of the fully re- 
laxed graphene sheet with the molecule adsorbed to it. 
The strength of the molecular doping is discussed in light 
of the density of states and the highest occupied and low- 
est unoccupied molecular orbitals (HOMO and LUMO) 
of the adsorbate. The position of these orbitals, visible as 
peaks in the DOS, is practically independent of the ori- 
entation and adsorption site of the molecule, so we only 
show the total DOS for one geometry per molecule. We 
can now distinguish two charge transfer mechanisms: i) 
a charge transfer can occur due to the relative position in 
the DOS of the HOMO and LUMO of the adsorbate. If 
the HOMO is above the Fermi level of pure graphene (the 
Dirac point), there is a charge transfer to graphene. If the 
LUMO is below the Dirac point, charge will transfer to 
the molecule, ii) The charge transfer between adsorbate 
and graphene is also partially determined by the mix- 
ing of the HOMO and LUMO with the graphene orbitals 
(hybridisation). This mixing scales with the overlap of 
the interacting orbitals and the inverse of their energy 
difference. 

It is more difficult to discuss the adsorption energy in 
this way because of the large amount of possible inter- 
acting orbitals present in graphene. Our investigation 
starts with the non-magnetic molecules H2O, NH3, and 
CO, followed by the paramagnetic ones, NO2 and NO. 
We discuss the molecules in the order of increasing com- 
plexity of their charge transfer mechanism. 



A. H2O on graphene 

We examine the following orientations of the H2O 
molecule with respect to the graphene surface: starting 
from the O atom the H-0 bonds pointing up (u), down 
(d) or parallel to the graphene surface (n). Another ori- 
entation (v) was suggested in a theoretical study, based 
on an empirical method, of the adsorption of H2O on 
graphite. This orientation has one 0-H bond parallel 
to the surface and the other one pointing to the surface. 
All properties were found to be almost invariant with re- 
spect to rotations around the axis perpendicular to the 
surface and through the oxygen atom, and therefore we 
will not discuss this orientation. The results of the cal- 
culations are given in table [ij 

TABLE I: H2O on graphene: the adsorption energy (^a), 
the distance of H2O above the graphene surface (c/), and the 
charge transfer from the molecule to graphene (AQ) for ten 
different geometries. 

Position Orientation ^a(nieV) d{k) AQ(e) 



B 


u 


18.4 


3.70 


0.021 


T 


u 


18.7 


3.70 


0.021 


C 


u 


20.3 


3.69 


0.021 


B 


n 


23.8 


3.55 


0.013 


T 


n 


23.7 


3.56 


0.015 


C 


n 


26.5 


3.55 


0.014 


B 


d 


17.8 


4.05 


-0.009 


T 


d 


18.5 


4.05 


-0.009 


C 


d 


19.4 


4.02 


-0.010 


C 


V 


47.0 


3.50 


-0.025 



From table |l] we learn that all the adsorption energies 
are small, which is partially a consequence of the calcu- 
lation method. This is not very important because the 
adsorption energies are only used to compare the differ- 
ent geometries and to find the best position and orien- 
tation of the molecule for which we need only relative 
values. Table 1 also shows that the adsorption energy 
is primarily determined by the orientation (u, d, n, v) 
and to a lesser degree by the position (C, B, T) of the 
molecule. The energy differences are 5-6 meV with re- 
spect to the orientation, but they vary with about 1-2 
meV when changing the position. This difference in im- 
portance of position and orientation is even more pro- 
nounced when we look at the charge transfers. If the O 
atom points to the graphene surface, there is a (small) 
charge transfer to graphene, but if the H atoms point to 
the surface there is a small charge transfer to the wa- 
ter molecule. This is a consequence of the form of the 
HOMO and LUMO of H2O and their relative position 
with respect to the Dirac point (see Fig. [2|. 

The HOMO (Ibi) is completely located on the O atom, 
but the LUMO (4ai) is mostly located on the H atoms. In 
the u orientation the HOMO plays the dominant role and 
donates, through a small mixing with graphene orbitals 
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FIG. 2: (Color online) H2O on graphene. Inset: (a) the 
HOMO and (b) the LUMO of H2O (the H atoms are white 
and the oxygen atom is red; green and yellow indicate differ- 
ent signs of the orbital wavefunction) . Main panel: DOS of 
H2O on graphene. The blue dotted lines show the positions 
of the molecular orbitals of H2O. 

above the Fermi level, some charge to graphene. There is 
also a (stronger) mixing with the orbitals below the Dirac 
point, because they are closer in energy, but this does not 
induce any charge transfer because all these orbitals are 
filled. In the d and v orientation, the LUMO of H2O 
interacts stronger with the surface and is able, through a 
small mixing with the graphene orbitals below the Dirac 
point, to accept some charge from graphene. There is 
also a stronger mixing with orbitals above the Dirac point 
now, but this does not induce any charge transfer because 
all these orbitals are empty. In the n orientation it is 
again the HOMO that will interact stronger, but now 
there is also some interaction with the LUMO. There 
will be a charge transfer from the molecule to graphene, 
but, because of the interaction with the LUMO, it will 
be smaller. 

Experimentall}E! one finds that H2O acts as an acceptor 
on graphene which is in accordance with our theoretical 
results where we find that the acceptor character (C,v) 
is energetically favored on perfect graphene. 



B. NH3 on graphene 

Two orientations of the ammonia molecule were inves- 
tigated, one with the H atoms pointing away from the 
surface (u) and the other with the H atoms pointing to 
the surface (d). All properties were again found to be 
almost invariant to rotations around the axis perpendic- 
ular to the surface and through the nitrogen atom. The 
results of the calculations are given in table [llj The ad- 
sorption site and the orientation are now seen to be of the 
same importance for the adsorption energy. The charge 
transfer, however, is solely determined by the orientation 
of the NH3 molecule. 

There is a small charge transfer from the molecule to 
the graphene surface of 0.03 e in the u orientation and 
there is (almost) no charge transfer in the d orienta- 
tion. We can see how this comes about by looking at 



TABLE II: NH3 on graphene: the adsorption energy (^a), 
the distance of NH3 above the graphene surface (d), and the 
charge transfer from the molecule to graphene (AQ) for six 
different geometries. 

Position Orientation ^a(meV) d{k) AQ(e) 



B 


u 


21.1 


3.86 


0.026 


T 


u 


20.1 


3.86 


0.026 


C 


u 


30.8 


3.81 


0.027 


B 


d 


14.7 


4.08 


0.001 


T 


d 


15.6 


3.97 


0.000 


C 


d 


24.7 


3.92 


-0.001 



the HOMO (3ai) and LUMO (4ai) of the NH3 molecule 
(Fig. |3|a) and (b)). In the u orientation, the HOMO is 
the only orbital that can have a significant overlap with 
the graphene orbitals and thus can cause charge transfer. 
As a consequence the NH3 molecule will act as a donor. 
In the d orientation, both HOMO and LUMO can cause 
charge transfers which are similar in magnitude but in 
opposite directions. The net charge transfer is therefore 
close to 0. 

The u orientation is energetically favored which explains 
the donor character as observed experimentall}^. 
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FIG. 3: (Color online) NH3 on graphene. Inset: (a) the 
HOMO and (b) the LUMO of NH3 (the N atom is blue and the 
H atoms are white). Main panel: DOS of NH3 on graphene. 

We also performed LDA calculations for the adsorp- 
tion of NH3 on graphene. The results of these are sim- 
ilar for the charge transfer, but the adsorption energy 
is much larger (±100 meV). The real adsorption energy 
is known to lie between the two approximate values ob- 
tained through GGA and LDA. This is consistent with 
the experimental observatiorPthat the adsorbates can be 
removed from the surface by annealing at 150° C. 



C. CO on graphene 

Three different orientations were used for the CO 
molecule. Two with the molecule perpendicular to the 
surface, with the O atom above the C atom (u) and the 
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other way around (d), and one parallel to the surface (n). 



D. NO2 on graphene 



TABLE III: CO on graphene: the adsorption energy (Ea), 
the distance of CO above the graphene surface (d), and the 
charge transfer from the molecule to graphene (AQ) for seven 
different geometries. 

Position Orientation ^a(nieV) d{A) AQ{e) 



B 
T 
C 
T 
C 
B 
C 



10.0 
9.6 

13.1 
8.4 
9.6 

14.0 

14.1 



3.75 
3.75 
3.73 
3.72 
3.70 
3.74 
3.74 



0.019 
0.019 
0.019 
0.009 
0.010 
0.013 
0.012 



From table III we notice that the CO molecule always 
acts as a donor. The size of the charge transfer only de- 
pends on the orientation of the molecule with respect to 
the surface. The differences in charge transfer are due to 
differences in orbital overlap between the HOMO (5cr) of 
the CO molecule and graphene. The LUMO (27r) seems 
to play no important role in the doping process although 
it is closer to the Dirac point than the HOMO. To under- 
stand this we have to take into account the symmetry of 
this orbital and the graphene orbitals. The DOS below 
and close (< 3 eV) to the Dirac point originates from 
(mostly) bonding combinations of the p^ atomic orbitals 
of the C atoms of graphene. The DOS above the Dirac 
point is mostly due to anti-bonding combinations. The 
completely anti-symmetric LUMO will therefore mostly 
interact with the DOS above the Dirac point which does 
not cause any doping. 

The HOMO is thus the more important orbital and the 
charge transfer is consequently always to graphene. Be- 
cause the HOMO is mainly located on the C atom, the 
charge transfer is the largest when the C atom is closest 
to the surface (u orientation), the smallest when the O 
atom is closer to the surface (d orientation ) and interme- 
diate when both atoms are at an (almost) equal distance 
from the surface (n orientation). 
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FIG. 4: (Color online) CO on graphene. Inset: (a) the HOMO 
and (b) the LUMO of CO (the C atom is black and the O atom 
is red). Main panel: DOS of CO on graphene. 



In Ref. [17] it was stated that adsorbates with a mag- 
netic moment in general result in a larger doping. There- 
fore we will now turn our attention to paramagnetic 
molecules. The first one is NO2 which has, in a spin- 
polarized calculation, an energy that is 0.4 eV smaller as 
compared to an unpolarized one and therefore is param- 
agnetic. We examine three different orientations of the 
NO2 molecule: starting from the N atom, the N-0 bonds 
pointing up (u), down (d) or parallel to the graphene 
surface (n). 



TABLE IV: NO2 on graphene: the adsorption energy (Ea), 
the distance of NO2 above the graphene surface (d), and the 
charge transfer from the molecule to graphene (AQ) for six 
different geometries. 



Position Orientation Ea (meV) 



B 


d 


67.4 


3.61 


-0.099 


T 


d 


65.3 


3.61 


-0.099 


C 


d 


62.6 


3.64 


-0.098 


B 


u 


54.7 


3.83 


-0.089 


T 


u 


54.5 


3.93 


-0.090 


C 


n 


66.7 


3.83 


-0.102 



The LUMO (6ai,i) of NO2 is located 0.3 eV below the 
Dirac point (Fig. [5| . This induces a large charge transfer 
to the molecule. But there are also some NO2 orbitals 
close enough to the Dirac point to cause some charge 
transfer in the opposite direction through orbital mixing 
(especially the HOMO (6ai,|)). The latter charge trans- 
fer is smaller than the first one but it is still noticeable in 
the strength of the magnetic moment of the system. The 
total magnetic moment of graphene and adsorbate in e.g. 
the (B, d) orientation is 0.862 /i^. The charge transfer 
from graphene (M=0/i5) to NO2 (M=l/i6) is 0.099 e, so 
the orbital mixing causes a charge transfer of ±0.039 e 
to graphene. 



E. NO on graphene 

To test whether or not there is always strong doping 
by paramagnetic molecules, we will investigate another 
one. For NO a spin polarized calculation gives an en- 
ergy that is 0.3 eV lower than a non-spin-polarized one, 
so NO is indeed a paramagnetic molecule. We investi- 
gate the same orientations and use the same notations 
as for the CO molecule (replace C with N). Contrary to 
the claim made in Ref. 17, we did not find that NO in- 
duces any strong doping. The charge transfers are an 
order of magnitude smaller than in the case of the NO2 
molecule (see table |V|) which is comparable to the non- 
magnetic molecules. Physically, we can understand this 
if we compare the DOS of the adsorbates NO (Fig. [6| 
and NO2 on graphene. For NO2 adsorbed on graphene. 
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FIG. 5: (Color online) NO2 on graphene. Inset: HOMO and 
LUMO of NO2 (the N atom is blue and the O atom is red). 
Main panel: spin-polarised DOS of NO2 on graphene. 



the LUMO is situated 0.3 eV below the Dirac point of 
graphene. This induces a strong doping. In the case of 
NO (Fig. [6| the HOMO is degenerate (27r^, 27r^) and 
is half filled (so it is also the LUMO) and lies only 0.1 
eV below the Dirac point. This induces a very small 
charge transfer from graphene to NO, but, due to its 
small strength, it can be (over) compensated by orbital 
mixing. The HOMO/LUMO of NO can, because it is 
half filled, cause charge transfer in both directions by 
mixing with the graphene orbitals below and above the 
Dirac point. But, as in the case of the LUMO of CO, 
it interacts mostly with the latter due to symmetry rea- 
sons. So the orbital mixing leads to charge transfer to 
graphene. 




Energy (eV) 



FIG. 6: (Color online) NO on graphene. Inset: (a) 5cr orbital 
and (b) HOMO/LUMO of NO (the N atom is blue and the 
O atom is red). Main panel: spin-polarized DOS of NO on 
graphene. 

and graphene leads to a larger orbital overlap and conse- 
quently to more orbital mixing (i.e. a larger charge trans- 
fer). The differences in the distance can be explained by 
the overlap of the 5<j orbital. The position of this or- 
bital is very close in energy to that part of the graphene 
DOS that originates from the (bonding) combinations of 
carbon p^ orbitals around the F point. Mixing of these 
orbitals induces a net energy shift upwards so they re- 
pel each other strongly. The geometry of the 5<j orbital 
gives a large overlap in the u orientation, a smaller over- 
lap in the d orientation and the smallest overlap in the 
n orientation. This gives a simple explanation for all the 
differences found from our calculations. 



TABLE V: NO on graphene: the adsorption energy (Ea), the 
distance of NO above the graphene surface (d) , and the charge 
transfer from the molecule to graphene (AQ) for six different 
geometries. 

Position Orientation ^a(meV) d{A) AQ(e) 



c 


u 


15.7 


4.35 


0.006 


T 


u 


14.0 


4.35 


0.006 


C 


d 


12.6 


4.11 


0.007 


T 


d 


10.6 


4.27 


0.005 


C 


n 


27.9 


3.71 


0.018 


B 


n 


28.5 


3.76 


0.017 



IV. SUMMARY AND CONCLUSIONS 

The charge transfer between the considered adsorbates 
and graphene is found to be almost independent on 
the adsorption site but it does depend strongly on the 
orientation of the adsorbate with respect to the graphene 
surface. We compared two paramagnetic molecules, 
NO2 and NO, and found that NO2 induces a relatively 
strong doping (-O.le), but NO does not (<0.02e). 
This is in contrast to Ref. [ITj where it was claimed 
that paramagnetic molecules are strong dopants which 
we found indeed to be the case for NO2 but not so for NO. 



We see from table [V| that the charge transfer due to 
orbital mixing always overcompensates the small trans- 
fer due to the position of the HOMO/LUMO, so NO 
always act as a donor. We notice that there are large 
differences in AQ and in the distance d. They obviously 
correlate because a smaller distance between adsorbate 



For the considered adsorbates the sign of the charge 
transf er a grees with what was found experimentally (see 
table VI) in Ref. [H Our results are also in good 



agreement with theoretical studies of the adsorption of 
molecules on large SWNT's in e.g. Ref. This suggests 
that some of the knowlegde of adsorption on nanotubes 
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TABLE VI: Summary of results. 



Adsorbate 


Theory 


Experiment^ 


Ea (meV) 


AQ(e) 


H2O 


acceptor 


acceptor 


47 


-0.025 


NH3 


donor 


donor 


31 


0.027 


CO 


donor 


donor 


14 


0.012 


NO2 


acceptor 


acceptor 


67 


-0.099 


NO 


donor 


/ 


29 


0.018 



should be transferable to graphene. 



Acknowledgments 

This work was supported by the Flemish Science Foun- 
dation (FWO-Vl), the NOI-BOF of the University of 
Antwerp and the Belgian Science Policy (lAP). 



* Electronic address: ortwin.leenaerts@ua. ac.bel 
^ Electronic address: bart.partoens@ua.ac.be| 

Electronic address: francois.p eeters@ua.ac.be| 
^ K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. 

Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov, 

Science 306, 666 (2004). 
^ K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. 

I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and A. A. 

Firsov, Nature (London) 438, 197 (2005). 
^ Y. Zhang, Y.-W. Tan, H. L. Stormer, and P. Kim, Nature 

(London) 438, 201 (2005). 
^ J. Kong, N. R. Franklin, C. Zhou, M. G. Chapline, S. Peng, 

K. Cho, and H. Dai, Science 287, 622 (2000). 
^ F. Schedin, A. K. Geim, S. V. Morozov, E. W. HiU, P. 

Blake, M. L Katsnelson, and K. S. Novoselov, Nature Ma- 
terials 6, 652 (2007). 
^ J. Zhao, A. Buldum, J. Han, and J. P. Lu, Nanotechnology 

13, 195 (2002). 

^ S. Santucci, S. Picozzi, F. Di Gregorio, L. Lozzi, C. Can- 

talini, L. Valentini, J. M. Kenny, and B. Delley, J. Chem. 

Phys. 119, 10904 (2003). 
^ J. A. Robinson, E.S. Snow, S.C. Badescu, T.L. Reinecke, 

and F.K. Perkins, Nano Letters 6, 1747 (2006). 
^ S. Peng, and K. Cho, Nanotechnology 11, 57 (2000). 
^° S. Peng en K. Cho, Nano Letters 3, 513 (2003). 

URL: http: //www. abini t.org/ [ 

J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77, 3865 (1996). 

B. S. Gonzalez, J. Hernandez- Roj as, J. Breton, and J. M. 
Gomez Llorente, arXiv:cond-mat/0706115 2vl (2007). 
N. Troullier, and J. L. Martins, Phys. Rev. B 43, 1993 
(1991). 

H. J. Monkhorst, and J. D. Pack, Phys. Rev. B 13, 5188 
(1971). 

F. L. Hirshfeld, Theor. Chim. Acta 44 , 129 (1977). 
T. O. Wehling, K. S. Novoselov, S. V. Morozov, E. E. 
Vdovin, M. L Katsnelson, A. K. Geim, and A. L Lichten- 
stein, arXiv:cond-mat/07Q33M| (2007). 



